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§ ■ ABSTRACT 

This paper analytically investigates a series of two-dimensional MHD reconnection 
solutions over a wide variation of magnetic Reynolds number (R* m ). A new series of 
solutions explains a continuous transition from Petschek-like fast regime to a Sweet- 
Parker-like slow regime. The inflow region is obtained from a Grad-Shafranov analysis 
used by Nitta et al. 2002 and the outflow region from a shock-tube approximation used 
by Nitta 2004, 2006. A single X-point (Petschek-like) solution forms for a sufficiently 
small R* em . As R* em gradually increases, the solutions shifts to an X-O-X solution with a 
magnetic island between two X-points. When R* m increases further, the island collapses 
to a new elongated current sheet with Y-points at both ends (Sweet-Parker-like). These 
reconnection structures expand self-similarly as time proceeds. As R* m increases, the 
reconnection rate and the reducible fraction of the initial magnetic energy of the system 
decrease as power-law functions of R* em - 
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1. INTRODUCTION electric resistivity. 

Majority believes that the magnetic reconnec- The entire scale ° f ^connection system is much 

i tion is an effective and universal energy conver- lar S er than typical particle scales, e.g., the ion 

sion process in many plasma systems. Especially, Larmor radius. This fact enables us to treat the 

KA i the magnetic reconnection plays essential roles in reconnection problems in a regime of magnetohy- 

*h" astrophysical phenomena, e.g., solar flares (see Adynamics (MHD). The MHD approximation is 

'• Tsuneta 1996), geomagnetospheric substorms, ex- v f ^ useful , to understand the basic macroscopic 

plosive phenomena in YSOs (see Koyama 1994, P^ 8108 , of tlie magnetic reconnection. Many the- 

Hayashi et al 1999) oretical and numerical works based on MHD sue- 

To apply the reconnection theory to astrophysi- ceeded to the , nature ° f the M ^ D recon- 

cal phenomena, the particular condition of the as- nectlon and Proposed its applications (see Sweet 

trophysical reconnection must be considered: It is 1958 > Parker c 1963 ' p etschek 1964, Vasyhunas 

intrinsically time- dependent and works in a "free- 1975 > Pri f 8t & F ° rbes 1986 ' a senes of P a P ers b ^ 

space" in which any environmental condition does U S ai and hl c s collaborators, e.g Ugai & Zheng 

not influence the evolution of the reconnection 2005 ' Sato & Ha y ashl 1979 f c). The magnetic 

system. In such condition, the self-similar evolu- reconnection properly takes place m resistive plas- 

tion model (see Nitta et al. 2001; hereafter paper m as. We should remember here that the electric 

1) of the reconnection will be relevant. In this resistivity cannot be reduced from MHD equa- 

paper, I discuss the dependence of the self-similar tlons ' but from some microscopic particle scale 

reconnection model on the wide variation of the P 1 ^ 8108 ' Thls has been the essential problem m 
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reconnection studies based on the framework of 
MHD. 

It is widely accepted that localized large elec- 
tric resistivity enhanced in the middle of the cur- 
rent sheet is required in order to establish a high 
power (fast) reconnection (Biskamp 1986, Scholer 
1989, Yokoyama & Shitaba 1994). Baty et al. 
(2006) claim a somewhat different result; main- 
taining the Petschek-like fast reconnection does 
not require significantly localized large resistivity. 
However, even in numerical simulation by Baty et 
al. (2006), it initially needs a localized (but not 
so large) resistivity to set up a Petschek-like fast 
reconnection. These results may suggest that it 
needs a localized resistivity to determine the re- 
connection point in the elongated current sheet 
system at least in the onset stage. 

Thus we may state that, in MHD reconnection, 
the so-called anomalous resistivity is required to 
establish the fast reconnection. Because the ori- 
gin of the anomalous resistivity is a microscopic 
process and outside the scope of MHD study, we 
must represent it by relevant model including pa- 
rameter (s). Unfortunately we neither have iden- 
tified the definite microscopic elementary process 
of the anomalous resistivity nor have understood 
its typical value yet (see Coppi & Friedland 1971, 
Ji et al. 1998, Shinohara et al. 2001). The 
value of enhanced resistivity may depend upon 
several factors. Recent particle simulations for 
the LHD instability (one of the probable candi- 
date for microscopic process of the resistivity en- 
hancement) revealed that the resistivity crucially 
depends upon the parameters: the current sheet 
thickness, the ratio of the electron plasma fre- 
quency to the electron gyro frequency, the guide 
magnetic field strength (see, Tanaka et al. 2006). 
This result inspires us to study the influence of 
the variation of the electric resistivity even in the 
same object, e.g., solar flares. 

Hence, in this paper, I discuss the MHD recon- 
nection in a free space for a very wide variation 
of the electric resistivity. The key parameter of 
this study is the magnetic Reynolds number R* em 
defined in section 3.1 in accordance with the mag- 
nitude of the electric resistivity. 

For the case of sufficiently small magnetic 
Reynolds number (R* em ~ 15), I discussed the 
plasma-/? dependence of the self-similar model in 
Nitta 2004; hereafter paper 3. As a result, I had 
clarified that the reconnection system is insensi- 
tive to the variation of plasma-/?. The reconnec- 



tion rate R* has rather universal value R* ~ 0.05 
for (3 < 1 and R* ~ 0.034 for (3 > 1 if the re- 
sistivity is sufficiently large (see figure 2 of paper 
3, more precise result is in figure 2 of Nitta 2006; 
hereafter paper 4). 

Paper 4 discussed the magnetic Reynolds num- 
ber dependence for the particular case that the 
reconnection point is fixed at a point in the 
middle of the current sheet. In this treatment, 
only the case with single X-point can take place 
(the Petschek-like solution). When the magnetic 
Reynolds number increases, a remarkable con- 
verging inflow concentrates into the region near 
the diffusion region. The direction of the con- 
verging inflow gradually tends to be parallel to 
the magnetic field lines as the magnetic Reynolds 
number increases. This converging inflow di- 
rectly results in the reduction of the reconnec- 
tion rate R* because it is proportional to the elec- 
tric field (the vector product of the velocity and 
the magnetic field) around the diffusion region. 
Thus we obtained the result R* oc -R* m _1 . How- 
ever, the converging inflow becomes too strong 
for R* em > 30: the converging speed exceeds the 
Alfven speed Vao at the asymptotic region (note 
Vao is the maximum limit for spontaneous flow 
speed). This is an irrational result and showing 
the limit of validity of the treatment in paper 4. 
The cause is likely from the approximation that 
the locus of the reconnection point is fixed. This 
deficiency must be improved to be applicable for 
the case of higher magnetic Reynolds number. 

The purpose of this paper is to clarify how large 
is the reconnection rate and how is the depen- 
dence on the reconnection structure at large mag- 
netic Reynolds number region in a more general 
circumstance. The hint is in the result of pa- 
per 4: the converging inflow in a large magnetic 
Reynolds number makes a channel flow structure 
(see section 4.1 of paper 4). From the discussion 
of the channel flow structure, we can speculate 
the following new type of solutions. 

In the channel flow, the inflow speed toward the 
slow shock is much larger than that toward the 
diffusion region. This non-uniform inflow leads a 
non-uniform compression of the field reversal re- 
gion. Since the compression outside the diffusion 
region is much stronger than that in the diffusion 
region, the reconnection point may move outside 
the present diffusion region. Thus, an O-point will 
form around the original diffusion region. Conse- 
quently, a typical X-O-X structure will form (see 
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section 4.6 of paper 4). 

In order to check the above speculative idea for 
very large magnetic Reynolds number, I extend 
the shock tube approximation used in papers 3 
and 4 to more general case with moving X-point. 
The alterations from the previous works are listed 
in section 3.2. 

This paper is organized as follows. The self- 
similar reconnection model is briefly summarized 
in section 2. The numerical scheme to solve the 
shock-tube approximated model is shown in sec- 
tion 3. The main result is listed in section 4. The 
detailed meaning of the obtained results and its 
applications are discussed in section 5. 

2. SELF-SIMILAR RECONNECTION MODEL 

We must note that many cases of the actual 
magnetic reconnections in astrophysical systems 
usually grow over a huge dynamic range in their 
spatial dimension. For example, the initial scale 
of the reconnection system can be defined by the 
initial current sheet thickness, but this is too 
small to be observed in typical solar flares. We 
do not have any convincing estimate of the scale, 
but if we estimate it to be of the order of the ion 
Larmor radius, it is extremely small (~ 10° [m] in 
the solar corona). Finally, the reconnection sys- 
tem develops to a scale of the order of the initial 
curvature radius of the magnetic field lines (~ 10 7 
[m] ~ 1.5% of the solar radius for typical solar 
flares). The dynamic range of the spatial scale is 
obviously huge (~ 10 7 for solar flares). For ge- 
omagnetospheric substorms, their dynamic range 
of growth is also large (~ 10 4 for substorms). 

In addition, we must note that, in many cases, 
the environmental timescale of energy storage 
(the convection timescale; typically of the order 
of day ~ 10 5 [sec] for solar active regions) is much 
longer than the evolution timescale of the recon- 
nection (the Alfven transit timescale; typically of 
the order of 10 2 [sec] for solar flares). This is com- 
mon with the geomagnetospheric substorms. 

Above mentioned properties (very wide dy- 
namic range of the expansion and the discrep- 
ancy of the timescales between the external cir- 
cumstance and the reconnection system) suggest 
that the external circumstance is almost station- 
ary during the evolution of the reconnection sys- 
tem. This means that the evolution of the mag- 
netic reconnection should be treated sponta- 
neous development in a "free space" , and that the 
external circumstance does not affect the evolu- 



tionary process of magnetic reconnection, at least 
at the expanding stage just after the onset of re- 
connection. We must study the spontaneous de- 
velopment of the magnetic reconnection for the 
astrophysical applications. 

Let us suppose the initial state to be a two- 
dimensional equilibrium state with anti-parallel 
magnetic field distribution, as in the Harris solu- 
tion. When non-uniform magnetic diffusion takes 
place in the current sheet by some localized re- 
sistivity, magnetic reconnection will start, and a 
pair of reconnection jets is ejected along the cur- 
rent sheet (see section 2 of paper 1). This causes 
a decrease in total pressure around the reconnec- 
tion point. Such information propagates outward 
as a rarefaction wave. In a low-/5 plasma (/? C 1 
in the region very distant from the current sheet 
[asymptotic region]; as typically encountered in 
astrophysical problems), the propagation speed 
of the fast-magnetosonic wave is isotropic, and 
is much faster than the slow mode wave. Thus, 
information about the decreasing total pressure 
propagates almost isotropically fast-mode 
rarefaction wave (hereafter FRW) with a speed al- 
most equal to the Alfven speed Vao in the asymp- 
totic region. Hence, the wave front of the FRW 
(hereafter FRWF) has a cylindrical shape except 
near the point where the FRWF intersects with 
the current sheet. The region swept by FRWF 
induces the inflow toward the reconnection re- 
gion owing to the total-pressure difference by the 
FRW. When the FRWF sufficiently expands, the 
initial thickness D of the current sheet becomes 
negligible comparing with the system size Vao^, 
where t is the time from the onset of reconnection 
(VAot 3> -Do)- 111 such a case, there is only one 
characteristic scale, i.e., the radius of the FRWF 
(VAot), which linearly increases as time proceeds. 
This is just the condition for self-similar growth. 

When the inflow toward the current sheet devel- 
ops sufficiently, a pair of slow shocks forms along 
the current sheet (the Petschek-like slow shock). 
After the formation of the slow shock, the en- 
ergy conversion drastically proceeds. The most 
of the magnetic energy is converted on the slow 
shock. Hence this stage is a "fast reconnection". 
Once this system of fast reconnection is set up, 
the dimension of the system unlimitedly develops 
self-similarly (see figures 5-7 of paper 1) until the 
FRWF reaches the external spatial scale. 

3. NUMERICAL SCHEME 
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3.1. Concept 

The basic concept of the shock-tube approxi- 
mation is the same with that of paper 4. As 
we mentioned in paper 3, the reconnection out- 
flow in our spontaneous model has a structure in- 
volving several discontinuities like a shock-tube 
problem (see figure 1 ). We solve the reconnec- 
tion outflow region by approximating as a kind of 
shock-tube problem (see details in the following 
subsections). In paper 4, we assumed a converg- 
ing inflow toward the slow shock, and discussed 
the structure of the outflow region as a function 
of the magnetic Reynolds number. Consequently, 
we obtained the magnetic Reynolds number de- 
pendence of the structure (see figure 7 of paper 
4) and the reconnection rate (see figure 6 of pa- 
per 4). 

In this paper we study the magnetic Reynolds 
number dependence of the outflow structure by 
the above shock-tube approximation in the regime 
of ideal (non-resistive) MHD in a more universal 
situation than that in paper 4. 

One may think that it is very curious that we 
can treat the magnetic Reynolds number in the 
regime of ideal MHD in this work. We should 
note, however, that if we can evaluate the re- 
connection electric field E z (note that E z is uni- 
form in the vicinity of the reconnection point 
from the Ampere's law because this region is al- 
most stationary in our self-similar model) and 
the magnetic field B* around the diffusion re- 
gion from the shock-tube problem, we can esti- 
mate the inflow speed v* (= E z /B*) toward the 
diffusion region. This inflow speed is in balance 
with the magnetic diffusion speed v* di ^ at the dif- 
fusion region. Providing that the Alfven speed 
Vaq in the asymptotic region to be fixed as a 
normalization value as in paper 3, inverse of the 
inflow speed toward the diffusion region shows 
the "effective" magnetic Reynolds number R* em 
{R* em = Vag/v^; = Vao/\v*\ , hereafter we abbre- 
viate "effective"). Thus we can treat the magnetic 
Reynolds number in a context of ideal MHD. 

Here we should note the difference of the defi- 
nition of the magnetic Reynolds number R* em = 
Vaq/v^j (where v* dif = rj/D, rj is the resistivity, 
D is the actual thickness of the diffusion region) 
from that of usual one (R em = Vao/[t]/L]) where 
L is the system scale. The value of R* em is consid- 
erably smaller than R em (R* em ~ 10~ 7 i? em ) since 
L ~ 10 7 D in typical solar corona. 



3.2. Difference from previous works 

The main frame of the logic is almost the same 
with our previous works (papers 3 and 4), but 
there is one essential difference: Supposing a free- 
dom of the motion of the reconnection point along 
the current sheet as follows. 

In our previous works, the reconnection point is 
presumably fixed at the center of symmetry which 
is in the middle of the current sheet. In this case, 
very strong converging inflow forms as R* em in- 
creases, and results in a channel flow structure 
(see figure 7 of paper 4). Finally, the horizontal 
component of the converging speed exceeds the 
asymptotic Alfven speed Vaq (see figure 4 of pa- 
per 4). This is cure because Vaq is the maximum 
speed of the spontaneous flow when all the initial 
magnetic energy converts to the kinetic energy. 
The reason of this curious result is probably in 
the above presumable assumption. 

In order to avoid this difficulty, let us introduce 
the freedom of the motion of the reconnection 
point (see section 1 and subsections 4.1 and 4.6 
in paper 4). This improvement enables us to nat- 
urally obtain a plausible result: the continuous 
transition of the structure as a result. 

3.3. Model 

We present a schematic picture of the reconnec- 
tion outflow in our self-similar evolution model in 
figure 1. First we introduce two coordinate sys- 
tems: the fixed coordinate and the zoom-out co- 
ordinate. The fixed coordinate is the conventional 
Cartesian coordinate in which a position vector r 
from the origin pointing a point at rest relative 
to the initial plasma in equilibrium is a constant 
vector. The zoom-out coordinate is defined as 
r' = r/[VAoi]- I n the zoom-out coordinate, any 
self-similarly expanding state looks to be station- 
ary. The coordinate axes of the zoom-out coor- 
dinate are defined as follows: x-axis is parallel 
to both the current sheet and initial anti-parallel 
magnetic field, y-axis is perpendicular to the cur- 
rent sheet, and z-axis is parallel (the right-handed 
direction) to the current sheet, but is perpendic- 
ular to the initial magnetic field (hence, d z =0 in 
this two-dimensional problem). 

We impose a freedom corresponding to the mo- 
tion of the reconnection point. The locus of the 
reconnection point (coincides with the diffusion 
region) x = x y is moving along the x— axis at 
the speed x v Vao relative to the initial plasma in 
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equilibrium because the system is self-similarly 
expanding. The symmetry with respect to x— 
and y— axes allows us to treat only the region 
x, y > 0. The outflow is composed of two differ- 
ent plasmas which have different origins. These 
two plasmas are touch in contact at the point 
x = x c (the contact discontinuity). A plasmoid 
forms around the contact discontinuity. The rear- 
half region (x < x c : the reconnection jet) is filled 
with the reconnected plasma coming from outside 
the current sheet (the inflow region). The front- 
half region (x > x c : the plasmoid) is filled with 
the original current sheet plasma. 

From the linear analysis of the inflow region 
in paper 2, we obtained the result that the in- 
flow velocity is parallel to the y— axis in the fixed 
coordinate (see equations [23] and [24] in paper 
2). In order to match the outflow region to the 
linearized inflow region, we impose a presumable 
assumption: x— component of the inflow veloc- 
ity vanishes in the fixed coordinate (see figure 2). 
This directly results in v xp = —x y at the recon- 
nection point in the zoom-out coordinate. We can 
clearly understand that, in the zoom-out coordi- 
nate, an apparent convergence into the diffusion 
region occurs owing to the reconnection point mo- 
tion along the current sheet (see figure 3). We can 
expect that this apparent converging flow results 
in a reduction of the reconnection rate as in the 
case of real converging inflow in our previous work 
(see section 3 of paper 4). 

When x y > 0, the reconnected magnetic flux is 
stored also in the region |x| < x y . This magnetic 
flux forms a magnetic island structure between 
two reconnection points at x — ±x y . I paid scant 
attention to this region because we cannot find 
any clear discontinuity in this region (from our 
preliminary numerical simulation: not published 
yet) and hence we cannot formulate shock- 
tube problem. For simplicity, we suppose a uni- 
form distribution of the stored vertical magnetic 
flux in the island. 

The entire outflow is surrounded by a slow shock 
which has a complicated 'crab-hand' shape (see 
Abe & Hoshino 2001). The Petschek-like slow 
shock (an oblique shock) is elongated from the dif- 
fusion region with a slight opening angle 9. When 
x y = 0, the Petschek slow shock has figure-X 
shape. If x y > 0, the reconnection point separates 
into two points (x = rkx y ). There is a reverse fast 
shock (an almost perpendicular shock) inside the 
reconnection jet (x — xf) in some cases. In front 



of the plasmoid, a forward figure- V-shaped slow 
shock (an oblique shock) forms. The opening an- 
gle and the locus (the crossing point with x— axis) 
are and x s , respectively. 

The entire structure including the several dis- 
continuities is analogous to the one dimensional 
shock-tube problem. I approximate this reconnec- 
tion outflow as a quasi one-dimensional problem 
in order to solve it analytically. Such an approxi- 
mation may be valid near the x— axis, because the 
system is symmetric with respect to the 
Hence we focus our attention on the quasi-one- 
dimensional problem along the inflow stream line 
and the reconnection outflow as follows. 

I assume the non-converging inflow seen in the 
fixed coordinate in order to keep the consistency 
with the inflow region. This is the critical differ- 
ence from the previous work (paper 4). I treat 
the figure L-shaped region x = x y and y = (a fi- 
nite hight up from the current sheet [finite 3> D 
where D is the initial current sheet thickness]). 
Each region between two neighboring discontinu- 
ities is approximated to be uniform. 

We note about the up-stream region p just 
above the Petschek-like slow shock. This region 
between the slow shock and the separatrix field 
line (the critical field line reaching the reconnec- 
tion point: see SFL in figure 3) can be also ap- 
proximated to be uniform. In this region, each re- 
connected field line has an almost straight shape 
and crosses the shock while each field line has hy- 
perbolic shape in the region above the separatrix 
field line. 

The region 3 between the contact discontinuity 
and the V-shaped slow shock looks to be non uni- 
form. Since we do not know how we can solve 
the two dimensional structure of this region an- 
alytically, I roughly approximate that the total 
pressure is uniform even in this region. 

The quantities denoting the initial uniform 
equilibrium at the asymptotic region are gas 
pressure Pq, mass density po an d magnetic field 
strength B . The plasma-/? value at the asymp- 
totic region is defined as (3 = P / (Bq/[2/j]) where 
/i is the magnetic permeability of vacuum. In the 
rest of this paper, we use the normalization of 
physical quantities as like in paper 2. We de- 
fine units for each dimension as follows: (unit of 
velocity) = Vao = B / y/JIpo (Alfven speed at the 
asymptotic region), (unit of the length) = VAot 
where t is the time from the onset of reconnec- 
tion, (unit of mass density)= p , (unit of magnetic 
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field) = B , (unit of pressure)= ■ po V| . 

Let us set the quasi-one-dimensional shock- 
tube-like problem as follows. The system has 23 
unknown quantities: P p , p p , v xp , v yp , B xp , B yp , 9, 
Pi, pi, Ui, B u x f , P 2 , p 2 , v 2 , B 2 , P 3 , Ps, v y3 , B x3 , 
B y 3, (ft and x s where P*, p*, v*, P* denote the 
pressure, density, velocity, magnetic field, respec- 
tively (Note x c = v x s = v 2 because no mass flux 
passes through the contact discontinuity). The 
suffices p, 1, 2 or 3 denote the region divided by 
the discontinuities (see figure 1). The suffices x or 
y denote the vector components. Quantities 9, Xf, 
4> and x s denote the inclinations of the Petschek- 
like slow shock, the locus of the fast shock, the 
inclination of the V-shaped slow shock and the 
locus of the V-shaped slow shock (crossing point 
with x— axis), respectively. 

These unknowns should be related to each other 
via conditions coming from the integrated form 
of conservation laws (i.e., the Rankine-Hugoniot 
[R-H] conditions) or other relations. The set of 
relations is listed in the next subsection. 

3.4. Basic Equations 

There is no essential difference of the basic equa- 
tions for the outflow structure between the equa- 
tions used in paper 4 except that we must include 

According to the above model of the reconnec- 
tion outflow, we obtain the following 22 equations 
for 23 unknown quantities (Detailed forms of each 
equation are listed in the appendix). I choose v yp 
as the variable controlled by hand (see section 
3.5). 

1. relations between pre-Petschek-like slow shock 
and asymptotic region 

(a) Frozen-in condition (eq. [Al]) 

(b) Polytropic relation (eq. [A2]) 

2. R-H jump conditions at Petschek-like slow 
shock 

(a) Pressure jump (eq. [A3]) 

(b) Density jump (eq. [A4]) 

(c) Velocity jump (parallel comp.) (eq. [A5]) 

(d) Velocity jump (perpendicular comp.) (eq. 
[A6]) 

(e) Magnetic field jump (parallel comp.) (eq. 
[A7]) 

(f) Magnetic field jump (perpendicular comp.) 
(eq. [A8]) 



3. R-H jump conditions at reverse fast shock 

(a) Pressure jump (eq. [A9]) 

(b) Density jump (eq. [A10]) 

(c) Velocity jump (eq. [All]) 

(d) Magnetic field jump (eq. [A12]) 

4. Magnetic flux conservation at reconnection 
point (eq. [A 13]) 

5. Force balance at contact discontinuity (eq. 
[A14]) 

6. R-H jump conditions at forward V-slow shock 

(a) Pressure jump (eq. [A15]) 

(b) Density jump (eq. [A16]) 

(c) Velocity jump (parallel comp.) (eq. [A17]) 

(d) Velocity jump (perpendicular comp.) (eq. 
[A18]) 

(e) Magnetic field jump (parallel comp.) (eq. 
[A19]) 

(f) Magnetic field jump (perpendicular comp.) 
(eq. [A20]) 

7. Boundary Condition at the tip of the outflow 
(eq. [A21]) 

8. Magnetic flux conservation all over the outflow 
(eq. [A22]) 

We can solve the most of these equations by 
hand, and by substituting the solutions into other 
equations, we can reduce equations. Finally, ten 
equations [Al], [A3], [A4], [A5], [A6], [A7], [All], 
[A14], [A21] and [A22] remain as complicated 
nonlinear coupled equations for ten unknowns p p , 
v xp , B xp , B yp , 9, Pi, pi, Xf, (p and x s with a con- 
trollable variable v yp . 

3.5. Numerical procedure 

We solve these coupled ten equations by an iter- 
ative method (the Newton- Raphson method). In 
order to find the well converged solution of the 
highly nonlinear coupled equations, a precise ini- 
tial guess of the unknowns is required. In general, 
the step to find an appropriate initial guess is the 
core of difficulty (this should be comparable to a 
"treasure hunting in ten-dimensional space" with 
an incomplete treasure map). Fortunately, this 
most difficult step had already been cleared in 
paper 4 for the case x y = 0. 

We can start from the solution of the case 
x y = for an arbitrary value of Our most in- 
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terest is focused on low-/? cases typically encoun- 
tered in astrophysical problems. We demonstrate 
the numerical procedure for the case (3q = 0.01 as 
the reference case because this is a typical value 
in the solar corona and geomagnetosphere. In pa- 
per 4, we have already had the solution in the 
spontaneous maximum inhalation case (x y = 0) 
for fa = 0.01. 

First, we start with the solution of the case 
Vyp = —Vi nS p(f3o = 0.01) = —0.06398 (where 
Vinsp[Po — 0.01] is the spontaneous maximum in- 
halation speed for f3 = 0.01) as the initial guess 
for the case v yp = — v insp (0.01) — Av y , where Av y 
is the decrement of v yp (see paper 3). Once we 
find a converged solution of the Newton-Raphson 
procedure for the case v yp = — v insp (0.01) — Av y , 
we treat it as the initial guess for the 
— Vins P (0. 01) — 2Av y , and we have successively ob- 
tained the series of solutions for different v yp . Of 
course, Av y should be a small enough value to 
keep good convergence of the Newton-Raphson 
method. From several trials, I carefully adopt the 
decrement: Ai^ = 10~ 8 . 

We obtain the result that, as \v yp \ increases, the 
locus x y of the reconnection point drastically in- 
creases (see figure 12: this means that the re- 
connection point is moving in x— direction at a 
speed x v Vao)- Simultaneously the strength of 
the reverse fast shock decreases. At a critical 
value \v yp \ = —0.063975, the reverse fast shock no 
longer forms (the density jump ratio and the pres- 
sure jump ratio reduce to unity). At this point 
we exchange the numerical code to another one 
for the case with no reverse fast shock forms (see 
the last part of appendix), and we can continue 
the calculation. Above procedure is valid for any 
value of flo- 

The above scan of v yp continues as long as we 
can obtain effectively converged solutions. I per- 
formed the calculation with 16 digit degree of 
complex numbers. I set the convergence goal 
of the Newton-Raphson method to be 10 digit 
degree of precision. In this case, R* em ~ 2100 
(where R* em = B xp / (v xp B yp - v yp B xp ): the nor- 
malized £?x.B-drift velocity) is almost the upper 
limit of the magnetic Reynolds number for effi- 
cient successive calculation. If we try further cal- 
culation, the Newton-Raphson method requires 
much smaller Av y and much more digit degree 
of numbers. It is so time consuming that we can- 
not efficiently proceed to much higher R* m regime. 
Fortunately, since we have already obtained suf- 



ficiently settled asymptotic behavior as shown in 
the following, I believe the further calculation is 
not worthy. 

3.6. Summary of procedure 

We here summarize our procedure. We solve 
the outflow region approximated as a quasi-one 
dimensional shock tube problem along the figure- 
L shaped line (x = x y and y — 0). Each region be- 
tween neighboring discontinuities is approximated 
to be uniform for simplicity. The system is de- 
scribed by two parameters: the plasma-/? f3 in 
the asymptotic region and the magnetic Reynolds 
number R* em . We fix (3 = 0.01 typical value 
for the solar corona. This problem has 23 un- 
known quantities: P p , p p , v xp , v yp , B xp , B yp , 9, 
Pi, pi, vi, Bi, x f , P 2 , p 2 , v 2 , B 2 , P 3 , P3, v y3 , B x3 , 
B y 3, and x s . We treat controllable vari- 

able handled by our hand. We can obtain other 22 
quantities by solving 22 coupled equations (Al)- 
(A22) denoting junction conditions over disconti- 
nuities. 

The magnetic Reynolds number (the key pa- 
rameter of this paper) is estimated by normal- 
ized i£x.B-drift velocity as R* m = B xp /(v xp B yp — 
v yp B xp ) by using obtained solution for unknowns. 
In the followings, the results are shown as func- 
tions of R* em . 

4. RESULT 
4.1. Reduction of the reconnection rate 

Our major interest is focused on the reconnec- 
tion rate. Figure 4 shows the magnetic Reynolds 
number {R* em ) dependence of the reconnection 
rate (R*) derived from the above approximated 
shock-tube problem in a very wide dynamic range 
of R* em (15 < R* em < 2100) where R* is de- 
fined as the normalized reconnection electric field 
R* = v xp B yp - v yp B xp in the unit of V A0 B . 
We can clearly find the asymptotic dependence 
R* oc (R* em )-\ 

We should note again here that the definition 
of the magnetic Reynolds number R* m adopted in 
this work is different from usual one: the quan- 
tity having dimension of the length included in the 
definition is not the system scale but the thickness 
of the diffusion region (see section 3.1). Hence the 
value of R* em is almost seven orders of magnitude 
smaller than usual one for the same resistivity if 
we estimate the thickness of the diffusion region 
to be of the order of the ion Larmor radius (~ 10° 
[m] in typical solar corona). 
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This reduction of the reconnection rate is not 
owing to the reduction of the inflow speed v yp , 
but to the enhancement of the apparent conver- 
gence of the inflow toward the reconnection point 
arising from the motion of the reconnection point 
(see figure 5). Even though the vertical inflow 
speed \v yp \ is almost constant, the apparent hori- 
zontal inflow speed (equals to the moving speed of 
the reconnection point) \v xp \ increases as R* m in- 
creases. The direction of the inflow velocity tends 
to be parallel to that of the magnetic field lines 
as R* em increases (see figure 6). This results in 
a reduction of the reconnection electric field, thus 
the reconnection rate reduces (remember that the 
reconnection rate is non-dimensional normalized 
value of the reconnection electric field). 

The reduction mechanism of the reconnection 
rate is almost the same with that of my previous 
work (see paper 4), but there is one crucially dif- 
ferent point. In paper 4, the actual converging 
speed |i> xp | exceeds the ambient Alfven speed Vao 
when R* em sufficiently increases. This is unrealis- 
tic because the spontaneous inflow speed cannot 
exceed the Alfven speed. The apparent converg- 
ing speed of this work asymptotically settles to 
the ambient Alfven speed for very large R* m - We 
must note that the direction of the magnetic field 
and that of the inflow velocity tends to be paral- 
lel even if \v xp \ is bound for the ambient Alfven 
speed (see figure 6). 

We find that the motion of the reconnection 
point results in the reduction of the reconnection 
rate. This means that the variation of the re- 
connection rate is owing to the transition of the 
structure of the reconnection system. This is the 
most important property of this model with mov- 
ing reconnection point. All of these changes are 
caused by the variation of the magnetic Reynolds 
number R* m . 

4.2. Continuous transition of the inflow 
structure 

This work clarifies the continuous transition of 
the entire structure from the Petschek-like one 
having large reconnection rate to the double re- 
connection point one having small reconnection 
rate in a series of solutions. When the magnetic 
Reynolds number is small enough {R* m ~ 15.63 
for /3 Q = 0.01: hereafter our discussion is limited 
to the representative case /3q = 0.01 because this 
is a typical value in the solar corona), a single 
reconnection point forms at the origin of the co- 



ordinate. When R* em > 15.63, the locus x y of the 
reconnection point shifts to the both directions 
along the outflow and the system forms double 
reconnection points at x — ±x y . Also a magnetic 
island forms between them (see figure 12). Note 
that we do not solve the region \x\ < x y by the 
shock-tube analysis. Instead of it, we assume a 
uniform distribution of the reconnected magnetic 
flux in this region for simplicity (see section 3.3 ). 

The shock-tube analysis provides the boundary 
condition to solve the inflow region by the Grad- 
Shafranov (G-S) method (see sections 3.2 and 3.3 
of paper 2). We impose this boundary condition 
on y = for simplicity. Of course, strictly to say- 
ing, this boundary condition should be imposed 
along the edge of the reconnection outflow, not 
on y = 0. The reasons to accept this simplicity 
are the followings: 1) The shape of the edge is 
rather complicated to be treated by simple way, 
2) Because the boundary condition obtained here 
is based on a drastic simplification (the quasi-one 
diminsional shock-tube approximation), it should 
be meaningless even if only here pursues strict- 
ness. 

By solving the G-S equation by a SOR routine 
(see paper 2), we can obtain the inflow struc- 
ture corresponding to individual outflow solu- 
tion obtained above. The resultant inflow struc- 
ture continuously changes according to the vari- 
ation of R* m : A typical Petschek-like structure 
{Rem ~ 15.63: see figure 7), a X-O-X type (figure 
8) or a double-Y point type solution (figure 9) for 
R* em > 15.63. We must note that only the region 
x 2 + y 2 < 1 is solved by the above G-S method 
in figures 7-9. The region x 2 + y 2 > 1 keeps its 
initial condition (an anti-parallel uniform equilib- 
rium configuration) and is not influenced by the 
reconnection yet. 

Whether the region \x\ < x y is occupied by an 
island or a collapsed current sheet-like structure 
is roughly determined by the ratio of the recon- 
nected magnetic field strength ([reconnected mag- 
netic flux]/ x y ) to the asymptotic magnetic field 
strength (-Bo)- If the ratio is greater than unity, 
an island forms, else (the ratio Cl) a collapsed 
current sheet-like structure forms. The critical 
magnetic Reynolds number is roughly estimated 
as R* em ~ 16.9 from the shock-tube calculation. 

4.3. Reducible fraction of the initial magnetic 

energy 



S. Nitta 



9 



The solutions for R* m > 15.63 include a mag- 
netic island or a current sheet-like structure (ex- 
tremely elongated and collapsed island) having a 
finite amount of magnetic energy. This means 
that the reconnection can release only a finite 
fraction of the total magnetic energy of the ini- 
tial current sheet system by one evolution of the 
self-similar reconnection. 

Let us estimate the reducible fraction of the ini- 
tial magnetic energy by single reconnection. Con- 
sider the separatrix stream line (SSL in figure 
10) reaching the reconnection point in the zoom- 
out coordinate (including the apparent converging 
flow). We approximate the shape of the separa- 
trix stream line to be straight for simple estima- 
tion. All the magnetic flux in the region below the 
separatrix stream line can reach the Petschek-like 
slow shock and its magnetic energy will be re- 
leased. We approximate here that the magnetic 
energy of the plasma which enters the Petschek- 
like slow shock is completely converted to the ki- 
netic / thermal energy of the reconnection outflow 
for simple estimation. By noting that this struc- 
ture is self-similarly expanding, we can calculate 
the reducible fraction of the initial magnetic en- 
ergy as follows. 

The system is initially just like a point at the 
origin in the conventional fixed coordinate. The 
system starts to expand self-similarly and sweeps 
the plasma in the initial current sheet system. 
At an instant, the region below a straight line 
(dashed line in figure 10) with inclination angle Or 
from the initial current sheet is the region already 
swept by the expanding past separatrix SSL or the 
region below the present SSL. Hence the magnetic 
energy in the triangle-shaped region (0 < < Or 
in figure 10) is reducible. We call this the "re- 
ducible region" . The angle Or is obtained as 

Or = tan _1 [(l - Xy)(v yp /v xp )] . (1) 

Thus the reducible fraction / is obtained as 

/ = R /^/2) . (2) 

The magnetic Reynolds number dependence of / 
is shown in figure 11. We can see that the re- 
ducible fraction drastically decreases as the mag- 
netic Reynolds number increases (asymptotically 

/ oc iCT 1/2 ). 

The intuitive meaning of the reducible fraction 
/ is as follows. Let us consider the flow struc- 
ture in the view point of the zoom-out coordi- 
nate. Note that the apparent converging inflow 



appears even in the region outside the FRWF in 
which there is no real flow in the fixed coordi- 
nate (see figure 10). All the stream lines in the 
region outside the FRWF radially converge to- 
ward the origin. We should note that only the 
plasma at the lower latitude region (< 9 R : the 
reducible region) from the (the original 

current sheet) reaches the region below the SSL 
when this plasma is swallowed to the expanding 
FRWF. This plasma can enter the Petschek-like 
slow shock within a finite time. On the contrary, 
the plasma at the higher latitude region (> Or) 
does not enter the Petschek-like slow shock and 
accumulates into the magnetic island or the new 
current sheet (see figures 8 and 9 ). Thus the mag- 
netic energy only in the reducible region (< Or) 
can be released by the self-similarly expanding re- 
connection system. 

5. SUMMARY AND DISCUSSION 

5.1. Summary 

We had investigated a magnetic reconnection in 
a free-space which is free from any external influ- 
ences in a series of papers (see papers 1-4). This 
paper is an extension to more universal situation. 
The new point of this study is introducing a free- 
dom of the motion of the reconnection point along 
the initial current sheet. I found the following ma- 
jor results. 

1) Reconnection rate 

The reconnection rate R* is a decreasing function 
of the magnetic Reynolds number R* em . This ten- 
dency is a natural result of the motion of the re- 
connection point. At the lowest value of the mag- 
netic Reynolds number (R* m = 15.63 for the case 
(3q = 0.01), the single reconnection point forms 
at the origin. As R* em increases, the locus of the 
reconnection point (\x\ — x y in the zoom-out co- 
ordinate) increases (see figure 12). This means 
that the reconnection point is moving along the 
original current sheet at a speed x v Vao- This mo- 
tion of the reconnection point naturally leads an 
apparent convergence of the inflow toward the re- 
connection region in the zoom-out (i.e., the recon- 
nection point-comoving) coordinate (see figure 3). 
Thus the reconnection rate (the normalized elec- 
tric field) reduces because the directions of the 
inflow velocity and the magnetic field tend to be 
parallel in the zoom-out frame (see figure 6). 

2) Continuous transition of solutions 

When x y = 0, the solution is very similar to the 
original Petschek model and having the largest re- 
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connection rate (typical fast reconnection). When 
x y > 0, two reconnection points form at x — ±x y . 
Between these reconnection points, a magnetic is- 
land or a new current sheet (exactly to saying, 
it is an extremely elongated and collapsed mag- 
netic island) forms. Hence, the system has a typ- 
ical X-O-X structure or a double Y-point struc- 
ture. This transition of the structure has close 
relation to the reduction of the reconnection rate 
as discussed above. As R* em increases, x y also in- 
creases. The structure continuously changes to 
the X-O-X type or the double Y-point type having 
smaller reconnection rate (slower reconnection). 
We can conclude that the structure directly cor- 
responds to the speed of the energy conversion. 
These structures expand self-similarly until some 
external circumstance influences the spontaneous 
evolution of the reconnection system. 

3) Reducible fraction of magnetic energy 
The reducible fraction / of the initial magnetic 
energy of the current sheet system is restricted 
in relation to the structures of the reconnection 
system in the above discussion. As x y increases 
(as R* em increases), / decreases (see figure 11). 
When the system has the X-O-X structure or the 
double Y-point structure, these structures remain 
the magnetic energy in the island or the newly 
formed current sheet-like structure. This means 
that a reconnection having these structure can re- 
lease only a finite fraction of the magnetic energy 
of the initial current sheet by one reconnection 
even though the current sheet magnetic energy is 
fully a free energy. 

In the following part, let us discuss the phys- 
ical meanings and the applications of the above 
results. 

5.2. Reconnection rate 

We obtained the magnetic Reynolds number 
dependence of the reconnection rate as R* oc 
(R* em )~ l . This behavior of the reconnection rate 
is closely associated with the motion of the re- 
connection point. When x y > 0, the reconnec- 
tion point is moving at the speed x v Vao in the 
direction parallel to the initial current sheet. In 
this work, the inflow direction is supposed to be 
perpendicular to the initial current sheet in the 
fixed coordinate to keep consistency with the re- 
sult of our linear analysis in the inflow region (see 
paper 2). In the zoom-out (reconnection point- 
co moving) coordinate, the inflow is oblique owing 
to the apparent horizontal component of the in- 



flow velocity (see figure 3). This apparent motion 
works as the converging inflow in paper 4 (see 
section 3 and figure 5). Hence we call this the 
apparent converging inflow. 

Figure 5 shows that absolute value of the hori- 
zontal component of the apparent converging in- 
flow does not exceed the ambient Alfven speed 
Vao (normalized as unity). This is a drastic 
change from the previous work (see figure 4 of pa- 
per 4). Needless to say, the ambient Alfven speed 
is the maximum value when all the magnetic en- 
ergy is converted to the plasma bulk kinetic en- 
ergy. Hence Vao is the maximum speed of any 
spontaneous flow. The current result (\v xp \ < 1 
even if R* em 3> 1) is plausible to be a realistic 
solution. 

Note that the curve of figure 4 terminates at 
R* em = 15.63 for the case (3 = 0.01. The author 
has the following somewhat speculative intuition 
about the dynamic behavior of the diffusion re- 
gion that is based on our experience of a lot of 
numerical simulations in paper 1. If we impose 
a very large electric resistivity in order to real- 
ize a situation with a small magnetic Reynolds 
number ( R* em < 15.63 ), the thickness of the dif- 
fusion region increases as time proceeds, thus the 
actual diffusion speed cannot increase almost at 
all. On the contrary, if we impose a very small 
resistivity to realize a situation with a large mag- 
netic Reynolds number, the thickness of the dif- 
fusion region decreases as time proceeds (if the 
mesh size is sufficiently small). The actual dif- 
fusion speed cannot decrease almost at all. This 
may show a self-regulation mechanism of the dif- 
fusion speed. The resultant reconnection rate is 
very insensitive to the variation of the imposed 
resistivity. Such dynamical response of the diffu- 
sion region had been empirically confirmed in our 
numerical MHD simulations. 

This dynamic behavior is mainly determined 
by the spontaneous inhalation speed by the 
Petschek-like slow shock (see section 5.4 of paper 
3). The maximum spontaneous inhalation speed 
is presented in figure 2 of paper 4. The value 
for Pq = 0.01 is \v yp \ = 0.06398 corresponding to 
R* em = l/\v yp \ = 15.63 (roughly to saying, this 
corresponds to 10 6 times the Spitzer resistivity or 
larger). This means that the magnetic diffusion 
speed cannot exceed 0.06398V^o- Even if the re- 
sistivity is very large, the diffusion speed will be 
self-regulated to this maximum value by sponta- 
neously adjusting the thickness of the diffusion 
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region. This clearly shows a passive nature of the 
diffusion region (see 4.2 of paper 4). 

This speculative discussion is consistent with 
the logarithmic (extremely weak) dependence of 
the reconnection rate on the magnetic Reynolds 
number in the original Petschek model. Such 
passive behavior of the diffusion region is ac- 
tual for the case of sufficiently large resistivity 
in which the diffusion region thickness is much 
larger than microscopic proper scales, e.g., the 
ion Larmor radius. In this case, the actual diffu- 
sion speed is almost constant independent of the 
resistivity, hence the magnetic Reynolds number 
R* em (= 15.63) (note that the definition of R* em is 
different from usual one; see section 3.1 ) and the 
reconnection rate (R* = 0.05) do not vary. This 
corresponds to the terminal point in the upper left 
end of the curve in figure 4. However, when the 
resistivity is very small and the diffusion thickness 
is comparable to its minimum value (e.g., the ion 
Larmor radius), the response of the reconnection 
system drastically changes to another scheme dis- 
cussed in section 4.1. 

On the contrary, if we force to enhance the in- 
flow speed v yp , we obtain the curious results with 
x y < for \vyp\ > 0.06398. This is nonsense be- 
cause it does not satisfy the realistic condition 
x y > 0. This may mean the possibility of the di- 
verging inflow, but this is also nonsense because 
the spontaneous inhalation cannot make any di- 
verging inflow by the fast-mode rarefaction (see 
Vasyliunas 1975). Thus I believe that we must 
terminate the curve in figure 4 at R* em (= 15.63) 
as like in this figure. 

We should note that such behavior in very large 
magnetic Reynolds number cases are very hard to 
be studied by MHD numerical simulations with 
finite sized mesh method. If we try to simulate 
this case, the diffusion region thickness thins to 
be almost equivalent to the mesh size. This is 
obviously nonsense because unrealistic numerical 
diffusion dominates the physical diffusion. This 
shows the limit of application of the MHD nu- 
merical simulation. 

The author is partially skeptical to the original 
Petschek's discussion on the logarithmic depen- 
dence of the reconnection rate to the magnetic 
Reynolds number in extremely high magnetic 
Reynolds number region. Petschek might pre- 
sumably assume that the diffusion region thick- 
ness can be infinitesimally thin in extremely high 
magnetic Reynolds number region to keep self- 



regulation process of the diffusion speed. We must 
note, however, that in the actual current sheet 
system, the diffusion region thickness is bound 
at the smallest value corresponding to the mi- 
croscopic proper scale (e.g., the ion Larmor ra- 
dius) . After the diffusion region thickness reaches 
its minimum value (of the order of the ion Lar- 
mor radius), the reconnection rate may not keep 
its value (as like in the very weak logarithmic de- 
pendence) but may decrease more rapidly (may 
be a power law) as like in the result of this work. 

5.3. Continuous transition from fast regime to 
slow regime 

We have obtained a wide variety of the solu- 
tions which have different values of x y correspond- 
ing to the variation of R* em . The transition of 
the structure results in the variation of the re- 
connection rate (see figure 4). The most signif- 
icant new point of this work is is to clarify the 
continuous transition from the fast regime to the 
slow regime. If the resistivity is large enough 
{R* em = 15.63), the Petschek-like fast regime takes 
place (reconnection rate R* ~ 0.05). As the resis- 
tivity decreases, the structure changes to the X- 
O-X type and finally reaches the double Y-point 
type (R* ~ 10~ 7 ). 

Another work which argue the continuous tran- 
sition from fast regime to slow regime is Priest 
& Forbes (1986). In their case, the reconnection 
solutions compose a family which is parametrized 
by constants denoting the external boundary con- 
dition at the inflow boundary (see section 3 of 
their paper). Hence their case is clearly based on 
the so-called "driven" model. On the other hand, 
in our case, the solution changes continuously ac- 
cording to the magnitude of the electric resistiv- 
ity which may relate to intrinsic properties of the 
plasma in the system. This work first clarifies the 
continuous transition based on the spontaneous 
model. 

Although the inverse scaling of the reconnec- 
tion rate R* with R* m is of the type normally as- 
sociated with super-slow reconnection (Priest & 
Forbes 2000), one must keep in mind that R* m is 
an effective magnetic Reynolds number based on 
the thickness of the diffusion region and not the 
global scale of the configuration. The reconnec- 
tion here is truly fast for small R* em because the 
electric field at the reconnection point is ~ 0.05 
of Alfven electric field defined by VaqBq, and is 
many orders of magnitude greater than would be 
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the case for super-slow or even slow reconnection. 
If we extrapolate our result to extremely large 
magnetic Reynolds number cases, R* decreases to 
~ 10~ 7 for the value of R* m corresponding to the 
Spitzer resistivity. This is just a typical value for 
the slow regime, i.e., the Sweet-Parker model. 

This drastic decrement of the reconnection rate 
is closely related to the transition of the struc- 
ture of reconnection system. Figures 7-9 clearly 
show the transition from the Petschek-like struc- 
ture (single X-point) to the double-Y structure 
via the X-O-X structure. During this transition, 
the locus of the reconnection points (x = ±x y ) 
separates in both direction along the initial cur- 
rent sheet at a finite speed (x y VAo)- The expected 
structure for extremely large magnetic Reynolds 
number is the double Y-point structure in which 
1) two reconnection points are moving in counter 
directions at the asymptotic Alfven speed V^o, 
and 2) a new current sheet-like structure forms 
between two reconnection points. Thus we can 
find a new solution for the slow reconnection 
regime which is morphologically very similar to 
the Sweet-Parker model especially in the inflow 
region. In addition, its reconnection rate (~ 1CT 7 ) 
is just similar to that of the Sweet-Parker model. 
We cannot distinguish them not only by the mor- 
phological structures but also by the energetics. 

However we should notice the critical differences 
between them. In our model, even in the case of 
small reconnection rate, 1) the diffusion region 
is localized around the Y-point (c.f., spread out 
into the entire current sheet in the Sweet-Parker 
model), 2) the magnetic energy is converted by 
the Petschek-like slow shock which elongates from 
the Y-point (c.f., converted by the magnetic dif- 
fusion in the entire current sheet in the Sweet- 
Parker model). When thinking from the view 
point of the reconnection rate, this is typical slow 
reconnection like the Sweet-Parker solution. How- 
ever, when thinking from the view point of the en- 
ergy conversion mechanism, it is considered as an 
extreme limit of the smooth transition from the 
fast regime as the Petschek solution, because the 
magnetic energy is mainly converted by the slow- 
mode wave (not the magnetic diffusion). There- 
fore, we must consider this is a quite new type of 
slow reconnection solution. 

5.4. Reducible fraction of magnetic energy by 
one reconnection 



The majority may think that the magnetic en- 
ergy stored in the current sheet system can be 
fully released by a reconnection. However, in 
section 4.3, we have clearly shown the magnetic 
Reynolds number dependence of the reducible 
fraction of the initial magnetic energy. If the 
electric resistivity decreases, released fraction of 
the magnetic energy by single reconnection also 
decreases (see figure 11). We can find that the 
magnetic energy of the current sheet system is 
analogous to the free energy of thermodynamic 
systems: the reducible fraction is restricted by 
additional conditions (e.g., the resistivity in this 
case) . 

Let us consider an ensemble of equivalent cur- 
rent sheet systems. While each configuration and 
stored magnetic energy is the same, the released 
energy can differ depending on the individual 
value of the resistivity enhanced in the diffusion 
region. If the resistivity is large enough and the 
Petschek-like reconnection takes place, the mag- 
netic energy will be completely released. On the 
contrary, if the resistivity is small and the double 
Y-point reconnection takes place, released energy 
will be very small. Finally the reducible fraction 
/ tends to / oc R* em ~ 1/2 . 

We should note that it is insufficient to estimate 
the magnitude of released energy by single recon- 
nection even if we can estimate the detailed con- 
figuration and the magnetic strength of the cur- 
rent sheet system. The released energy strongly 
depends on the magnetic Reynolds number which 
is very difficult to forecast beforehand. 

5.5. Observational estimation of reconnection 

rate 

In a lot of works, the reconnection rate is de- 
fined as the inflow Alfven Mach number. This 
definition, however, is not appropriate for the case 
having converging inflow as like in this study. We 
must note that the inflow Alfven Mach number 
itself does not denote the reconnection rate if the 
inflow is inclined with respect to the magnetic 
field lines in the zoom-out (reconnection point- 
co moving) frame. As an extreme example, if the 
inflow velocity is just parallel to the magnetic field 
lines, the reconnection rate is vanished even if the 
inflow speed is not vanished. Since the .ExiJ-drift 
determines the shift speed of the magnetic field 
line, the actual reconnection rate should be de- 
fined by the reconnection electric field not by the 
inflow speed: R* = v xp B yp — v yp B xp . The popular 
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estimation of the inflow Alfven Mach number for 
the reconnection rate coincides with true recon- 
nection rate only if v xp = (no converging flow 
or fixed reconnection point) and B xp = 1 (same 
with ambient magnetic field strength). 

Recent observation enables us to estimate the 
reconnection rate for actual reconnection events. 
Tsuneta (1996) first estimated the reconnection 
rate (the inflow Alfven Mach number ~ 0.07) in 
a solar flare event on 1992 Feb. 21 by Yohkoh 
observation. Yokoyama et al. (2001) estimated 
the inflow Alfven Mach number (~ 0.001 — 0.03) 
in a flare on 1999 March 18 by SOHO obser- 
vation. Isobe et al. (2005) and Nagashima & 
Yokoyama (2006) also estimated the inflow Alfven 
Mach number for several flares. If we note signif- 
icant ambiguity to estimate the Alfven speed, we 
may conclude that these estimated values of the 
reconnection rate are well consistent with almost 
theoretical upper-limit value (of the order of 10~ 2 ) 
of the fast reconnection models including the self- 
similar evolutionary model. 

However we should note that these studies can 
make sense to estimate the reconnection rate only 
if the inflow is approximately perpendicular to the 
magnetic field lines (realistically saying, it is plau- 
sible only if the reconnection point is at rest in 
the initial plasma- rest frame). If the reconnec- 
tion point is moving at a considerable speed com- 
paring with the Alfven speed, we must adopt the 
correct definition R* = v xp B yp — v yp B xp for esti- 
mation instead of the inflow Alfven Mach number. 
For example, this is important for the case of the 
loop-top flare events in which the reconnection 
point will rise upward at a considerable speed. 

5.6. Relation to particle acceleration 

The reconnection system may play an impor- 
tant role also for the particle acceleration. First, 
Tsuneta and Naito (1998) pointed out that the 
reconnection system having reverse fast shock in 
the outflow surrounded by the Petschek-like slow 
shock can work as an efficient electron accelera- 
tor. In their result, non-thermal electrons (20- 
100 KeV) can be produced in a very short time 
(0.3-0.6 sec) by the first order Fermi acceleration 
process around a quasi perpendicular fast reverse 
shock. A theoretical extension by Mann et al. 
(see Mann et al. 2006) argues the electron accel- 
eration in the same situation by the shock drift 
acceleration up to 10 MeV. 

We must note that, in the self-similar evolu- 



tionary reconnection model, the parameter range 
in which the reverse fast shock forms in the recon- 
nection outflow is not so wide (low [3: [3 < 0.441 
[see figure 2 of paper 4], not so high magnetic 
Reynolds number: R* m < 22.9 [see figure 4]). If 
the parameters are in this range, the self-similar 
reconnection system also can work as an effi- 
cient particle accelerator. Because the accelera- 
tion timescale (~ 10 _1 [sec]) estimated by Tsuneta 
& Naito (1998) is much shorter than the evolution 
timescale (the Alfven transit time ~ 10 2 [sec] for 
typical solar flares) of the reconnection system, 
the energy spectrum of the accelerated particles 
may be conserved during the evolution while the 
total amount (particle flux) of the accelerated par- 
ticles increases as the system scale increases in 
proportion to the time from the onset. 

5.7. Relation to micro/nano flares 

The property discussed in section 5.4 for the 
cases of very large magnetic Reynolds number 
(R* em 3> 1) inspires us a relation to micro- or 
nano-flares. As a result in section 4.3, the re- 
connection with very small resistivity can release 
only a small fraction of the magnetic energy of 
the initial current sheet system at one evolution 
even if the stored magnetic energy is huge. This 
may correspond to small energy-scale flares. 

As discussed in section 5.3, the release of the 
magnetic energy takes place on the slow shock 
elongated from the reconnection point. The ener- 
getic region (here after, the "hot region" ) around 
x ~ ±x y is moving along the initial current sheet 
at the speed ~ x v Vao where the locus x y (< 1) 
of the reconnection point is an increasing func- 
tion of R* em . These hot regions will be observed 
as a couple of two bright points (in X-ray or UV 
depending on the temperature) having a bipolar 
motion along the magnetic field. It will be seen 
as phenomena that bright points spread along the 
magnetic flux tube (typically seen in the so-called 
"loop flares", e.g., in the TRACE movies of C6.3 
flare on 2000 Sep. 5 or C6.1 flare on 2001 Aug. 
26). 

If the magnetic Reynolds number is somewhat 
larger (R* em ~ 15.6 — 20; see section 4.2) than 
that of the single X-point case, the remarkable X- 
O-X type reconnection occurs (see figure 8). The 
magnetic energy remained in the magnetic island 
cannot be released by further magnetic reconnec- 
tions because the remained island will be stable 
for magnetic reconnections. Hence this case is 
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one-off phenomena which can release a finite frac- 
tion of the magnetic energy (> 36% of the initial 
magnetic energy of the current sheet; from figure 
11). The restriction of the released fraction of the 
magnetic energy is related to the drastic topol- 
ogy change of the magnetic field lines (i.e., from 
"current sheet" to "island"). 

The very large magnetic Reynolds number cases 
{R* m 3> 20) which form the double Y-point type 
structure are more interesting. An extremely 
elongated magnetic island remains between the 
two Y-points. This collapsed island looks like a 
new current sheet which will be unstable again for 
reconnections; see figure 9. The author empha- 
size that this double Y-point type reconnection 
with moving hot regions can recurrently occur in 
a current sheet system. Even if the amount of 
released energy by one reconnection is extremely 
small, such recurrent occurrence of the double Y- 
point type reconnections finally can release a large 
fraction of the magnetic energy of the initial cur- 
rent sheet system. In addition, this type of re- 
connections may occur simultaneously to form a 
nested spatial self-similar structure. 

The first-light image of HINODE (Solar-B satel- 
lite) revealed a lot of very small X-ray bright 
points simultaneously take place on the entire 
solar surface. This phenomenon might corre- 
spond to the reconnections at very high magnetic 
Reynolds number. Such small, short time and 



faint events which we were not able to detect 
before HINODE might contribute to the coronal 
heating. In that case, we might be able to ob- 
tain the answer to the long term open question 
"coronal heating problem" soon. 
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APPENDIX 

EQUATIONS FOR STRUCTURE OF THE RECONNECTION OUTFLOW 
The structure of the reconnection outflow is determined by the following equations. 

We assume that the region between the asymptotic region and the pre-shock region is filled with non- 
resistive plasma. Hence, the magnetic flux is frozen into the induced inflow. We also assume a polytropic 
variation in the induced inflow because there is no violent process in the fast-mode rarefaction. Each 
region between neighboring discontinuities is approximated to be uniform. 

The essential difference from our previous works (papers 3 and 4) is introducing a new unknown x y 
(the locus of the reconnection point). However, we do not need new physical equation by assuming 
x y = —v xp (see section 3.4) in order to keep consistency with the linear perturbation theory in the 
inflow region. The alterations are in equations [A5], [A6], [A13], [A21] and [A22] as followings. 

Thus, we impose the frozen-in condition 

Bp _ B xp 
Po P P 

and the polytropic relation 

^OPO 7 = P pPp 



(Al) 



(A2) 
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where 7 is the specific heat ratio. We assume 7 = 5/3 (monoatomic ideal gas). 

There are several discontinuities in the reconnection outflow, i.e., X-shaped slow shock, reverse fast 
shock, contact discontinuity, and forward V-shaped slow shock (see figure 1). We set jump conditions 
for both sides of each discontinuity: First we consider the X-shaped slow shock Rankine-Hugoniot (R-H) 
jump conditions. The pressure-jump is 

p 

— — 1 + ■$-{— sin 9v xp + cos 9v yp ) 2 (X — 1) x - (cos 9B xp + sin 9B yp ) 2 /2 

x [-2Vj px X + (- sin 9v xp + cos9v yp ) 2 (X + l)]/[((- sm9v xp + cos9v yp ) 2 - XVl px f m }} ,(A3) 

where 



C-sp 



\ll P vlPp 



V Apx = yj (- sin 9B xp + cos 9B yp ) 2 / (pp p ) , 
p is the magnetic permeability of vacuum and X is the compression ratio. 



Then the density-jump is 



The velocity (parallel) jump is 



P — = X . (A4) 
Pp 



cos 0(vi - Xy) - v _ (— sin 9v xp + cos 9v yp ) 2 - V\ px 
sin 9v yp -v (- sin 9v xp + cos 9v 2 p - XVj px ) 



(A5) 



where v = (y yp B xp — B yp v xp )/ (B xp smO — B yp cos9) is the shift speed of the de Hoffmann- Teller 
coordinate. 
The velocity (perpendicular) jump is 

— sin 9{v\ — x y ) 1 



— sin 9v xp + cos 9y yp X 
The magnetic field (parallel) jump is 

sin 9Bi [(- sin 9v xp + cos 9v yp ) 2 - V\ px \X 



(A6) 



(A7) 



cos 9B xp + sin 9B yp (— sin 9v xp + cos 9v yp ) 2 — XV\ 
The magnetic field (perpendicular) jump is 

cos9Bi 
— sin 9B xp + cos 9B yp 

The compression ratio X is defined by the following equation (3rd order algebraic equation for X), 



{(— sm9B xp + cos 9 B yp ) 2 '[{cos 9 B xp + sm9B yp ) 2 {j — l)p p (— sm9v xp + cos9v. 



\2 



HI' I 



+ (- sin 9 B yp + cos 9 B yp ) 2 (27 P p - p p (— sm6v xp + cos9v yp ) 2 + 7P P (- sin 9v xp + cos 9v yp ) 2 )]}X 3 
+{— p p (— sin 9v xp + cos 9v yp ) 2 [(— sin 9B xp + cos9B y p) A {j — 1) 
+ {cos9B xp + sm.9B y p) 2 {j — 2)pp p (— sin9v xp + cos9v yp ) 2 
+ (— sin9B xp + cos 9 B yp ) 2 ((cos 9 B xp + sm9B yp ) 2 ('y + 1) 
+2p(2 r yPp — Pp(— sin 9v xp + cos9)v yp ) 2 + 7P P (— sin 9v xp + cos9v yp ) 2 )}}X 2 
+{pPp{— sm.9v xp + cos 9v yp [(cos 9 B xp + sm9B yp ) 2 'y + 2(— sin 9B xp + cos9B yp ) 2 (^ + 1) + p 
(2^P P - p p (— sin9v xp + cos 9v yp ) 2 + 7P P (- sin 9v xp + cos 9v yp ) 2 )]}X 
+{-(7 + 1 )p 2 Pp(- sin 9v xp + cos 9v yp f] = . 
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This equation has three roots. We must choose a real positive root larger than unity. 
The reverse-fast shock R-H conditions are as follows. The pressure-jump is 



(A9) 



where 
with 



C/ = 7 M 1 2 / (i-i/e / ) + (i-e?)M + i 

= H + \jl 2 + 2/A(2 - 7 )(7 + l) 7 M 1 2 / )/(2//? 1 (2 - 7)) , 
/ = 7 (1/A + 1) + (7 - 1)7^/2 , 
Pi = Pi/{Bl/{2^)) , 
M lf = (v 1 -XfV A )/c f , 



and 

The density-jump is 
The velocity jump is 

The magnetic field jump is 



c i = \/rPi/pi 

V A = B / y /JIp ~ . 



p2 -p 



v 2 -x f 



B 2 



(A10) 
(All) 

(A12) 



We must impose local magnetic flux conservation on both sides of region p and region 1. The magnetic 
flux conservation at X-point gives 



-V X pBy p + VypB Xp + (Vl - Xy)Bl = . 

Force balance at contact discontinuity is 

B\ B 2 x3 + 5 2 3 

P 2 + — = P 3 + — y - . 

2/i 3 2/x 

The V-shaped forward slow shock R-H conditions are as follows. The pressure-jump is 



_7_ 

c s0 



{x s smct)f{X h - 1) x 



J_ (B cos 0) 2 -2Vl 0x X h + (x s sin 0) 2 (X fe + 1) 



(A13) 



(A14) 



(A15) 



2 [(a;, sin 0) 2 - X h V2 0x ] 2 fip 

where c s0 = \fjPo~Jpo, Vaox — \J ( — sm 0A)) 2 / (mPo) an d ^ is the compression ratio. The density-jump 

P - = X h . (A16) 
Po 

The velocity (parallel) jump gives 



is 



cos <j)(v x3 - x s V A q) + sin 
-Vaqx s COS0 



7> 2 - T/ 2 

^mOx v AOx 



2 ' 



^mCte ~ ^hVAOx 



(A17) 
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where v m o x = x s sm<f) and Vao = -Bo/v^Po- The velocity (perpendicular) jump gives 

- sin <j)(v x3 - x s V A q) + cos <pv y3 _ 1 
V A qx s sin <f) X h 

The magnetic field (parallel) jump is 

cos <f)B x3 + sin <f)B y3 _ (v 2 mQx - V% 0x )X h 



(A18) 



(A19) 



5 O COS0 VmOx- X hV AQx 

The magnetic field (perpendicular) jump is 

-sin 0^,3 + cos (j)B y3 = i 2Q 
-B sin ' 1 ' 

The compression ratio X h is a solution of the following third order algebraic equation: 

{(B sin (p) 2 [(B sin 0) 2 (7 - l)p (^ sin 0) 2 + (5 sin (p) 2 (2-fP - p (x s sin 0) 2 + ~fp (x s sin 0) 2 )]}A^ 
+{-Po(^ sin0) 2 [( J B o sin0) 4 (7 + 1) + ( J B o cos0) 2 (7 - 2)pp (x s sin0) 2 
+ ( J B o sin0) 2 (( J B o cos0) 2 (7 + 1) + (2/i(2 7 P - p (x s sin0) 2 + 1Po (x s sin0) 2 )))]}X 2 
+{{ip 2 (x s sin 4>) 4 ((S cos0) 2 7 + 2(^0 sin 0) 2 (7 + 1) + i2(2<yP - p (x s sin 0) 2 + -fp (x s sin(f)) 2 ))}X h 

+{-( 7 +l)/i 2 p3( a : s sin0) 6 } = O. 

We must choose a real positive root larger than unity. 

The tip of the reconnection outflow touches the FRWF (see figure 1), hence A[ = at this point. 
This condition reduces to the following boundary condition at the tip of the outflow: 

—5^(1 — x c ) + -^[(l — x s ) tan0 — (x c — x y ) tan#] — B (l — x s ) tan0 = (A21) 

From magnetic flux conservation, the injected magnetic flux must be redistributed in the reconnection 
jet. This leads to the following equation for magnetic flux conservation: 

-B yp v xp + v yp B xp +[(x f - x y )B 1 + (x c - x f) B 2 \ = (A22) 
We assume the following trivial relations: 



Definition from the contact discontinuity 

x c = v 2 = v x3 

We can solve [A2] for P p , [A8] for B u [A9] for P 2 , [A10] for p 2 , [A12] for B 2 , [A13] for v yp , [A15] for P 3 , 
[A16] for p 3 , [A17] and [A18] for v x3 and v y3 , [A 19] and [A20] for B x3 and B y3 by hand, then substitute 
them into other nine equations ([Al], [A3], [A4], [A5], [A6], [A7], [All], [A14], [A21] and [A22]) for the 
following unknowns: p p , v xp , B xp , B yp , 6, P±, P i, Xf, <f) and x s . The only parameter included in this 
problem is the plasma-/3 value at the asymptotic region. By using a Newton-Raphson routine, with an 
initial guess for these unknowns, we obtain converged solutions. The procedure to obtain the series of 
converged solutions is discussed in section 3.4 in detail. 

As \vyp\ decreases, the strength of reverse fast shock reduces, and then the pressure jump Q and 
density jump £j simultaneously become unity at the critical value of \v yp \. At that point, the reverse 
fast shock vanishes. We convert the coupled equations to@another set compatible to the situation with 
no fast shock, i.e., several equations reform to the following equations: 
[A9] is replaced by 



P2 = Pl, 
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[A10] is replaced by 



P2 = Pi , 

[All] is replaced by 

vi=v 2 , 

[A12] is replaced to by 

B\ — B 2 , 

and [A22] becomes to be equivalent to [A13] and removed from the set of coupled equations. The 
number of coupled equations reduces to 21. This is consistent with that the locus Xf of the fast shock 
is no longer included in the set of unknowns and total number of unknowns becomes 21. 

APPENDIX 
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Fig. 1. — Schematic figure of the reconnection outflow. Several discontinuities form in the outflow. We consider a 
quasi-one dimensional shock-tube problem along the x— axis and x = x y . 

Fig. 2. — Schematic picture of the inflow structure observed in the "fixed" (initial plasma-comoving) coordinate. We 
assume the vertical inflow in order to keep consistency with our linearized analysis via Grad-Shafranov approach (Nitta 
et al. 2002). Note that the Alfven Mach number of this inflow does NOT mean the reconnection rate. 

Fig. 3. — Schematic picture of the inflow structure observed in the "Zoom-out" (reconnection point-comoving) 
coordinate. An apparent convergence (the horizontal component) of the inflow takes place. This horizontal velocity 
component increases as the magnetic Reynolds number R* em increases. The resultant oblique inflow tends to be parallel 
to the magnetic field lines as R* em increases. Thus the reconnection rate decreases as R* m increases. 

Fig. 4. — Magnetic Reynolds number dependence of the reconnection rate R* . This is the case for /3 = 0.01. The 
reconnection rate decreases as oc 1/R* m in very large magnetic Reynolds number. 
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Fig. 5. — Velocity and magnetic field at the region p versus magnetic Reynolds number R* em in the zoom-out 
(reconnection point-comoving) frame. This is the case for /?o = 0.01. These quantities are essential to determine the 
reconnection rate R* . As R* em increases, first, the transverse component B xp of the magnetic field quickly increases. 
Then, v xp drastically developed. The negative value of v xp denotes an apparent converging inflow toward the Petschek-like 
slow shock. Note \v xp \ does not exceed unity even in the range of very large magnetic Reynolds number. 

Fig. 6. — Direction of velocity and magnetic field at the region p in the zoom-out (reconnection point-comoving) 
frame. This is the case for (3q = 0.01. The circular dot and square denote v yp /v xp (tangent of the inclination angle of 
the inflow velocity) and B yp /B xp (tangent of the inclination angle of the inflow magnetic field), respectively. This graph 
clearly shows that the inflow velocity v p and the inflow magnetic field B p tend to be parallel to each other when the 
magnetic Reynolds number R* m increases. This causes the decrement of the reconnection rate R* = v p x B p . We should 
note that the information of the inflow AlfVen Mach number (v yp ) perpendicular to the current sheet is not sufficient to 
determine the reconnection rate. 

Fig. 7. — The magnetic field line configuration (Petschek type) in the inflow region for the case R* em — 15.63 and 
j3{) = 0.01 derived from linearized G-S approach. The single reconnection point forms at x — 0. The entire structure is 
very similar to that of the Petschek model. We call this the single X-point solution. The reconnection rate in this case 
is the maximum value R* = 0.0504. Note that only the region inside the FRWF (x 2 + y 2 < 1) is solved in this figure. 
The region outside the FRWF (x 2 + y 2 > 1) is in the initial uniform equilibrium. 

Fig. 8. — The magnetic field line configuration (X-O-X type) in the inflow region for the case R* em — 16.88 and 
/?o = 0.01. The double reconnection points are formed at x = ±0.048. A magnetic island forms between them. This is 
the X-O-X type solution. The reconnection rate in this case is R* = 0.0489. Note that only the region inside the FRWF 
(x 2 + y 2 < 1) is solved in this figure. The region outside the FRWF (x 2 + y 2 > 1) is in the initial uniform equilibrium. 

Fig. 9. — The magnetic field line configuration (double Y-point type) in the inflow region for the case R* m — 139.8 
and Pq — 0.01. A magnetic island forms between the double reconnection points at x = ±0.75. The island is extremely 
elongated and looks like a current sheet. This is a new type of slow regime. We call it the double Y-point solution. The 
reconnection rate in this case is R* = 0.00713. Note that only the region inside the FRWF (x 2 +y 2 < 1) is solved in this 
figure. The region outside the FRWF (x 2 + y 2 > 1) is in the initial uniform equilibrium. Although this solution looks 
like the Sweet-Parker solution, there are critical differences between them (see section 5.3). 

Fig. 10. — The schematic picture of the reducible region. The plasma in the region below the separatrix stream line 
(SSL) can enter the Petschek-like slow shock (SS) elongated from the reconnection point. Such plasma can release its 
magnetic energy. By noting the self-similar expansion of the system, the magnetic energy in the region below the dashed 
line having the inclination 8ji from the x-axis is reducible by single reconnection. We approximate the separatrix stream 
line to be straight and that the slow shock completely converts the magnetic energy for simplicity in order to estimate 
the reducible fraction. 

Fig. 11. — The reducible fraction of the magnetic energy of the initial current sheet system is shown as a function 
of the magnetic Reynolds number under the approximation that the separatrix stream line is straight. This is the case 
for p = 0.01. When the magnetic Reynolds number is small enough (R* m = 15.63), the reconnection can completely 
convert the magnetic energy (/ = 1). As the magnetic Reynolds number increases, the reducible fraction drastically 

i/2 

decreases. Finally it tends to / oc R* m 1 . This means that the released energy cannot be determined even if we 
estimate the magnetic energy of the current sheet system until we know the magnetic Reynolds number. 

Fig. 12. — The locus of the reconnection point (x = x y ) as a function of the magnetic Reynolds number R* em . This 
is the case for /3 = 0.01. We assumed vertical inflow in the fixed coordinate (thus x y = —v xp ). When the resistivity is 
large enough (the magnetic Reynolds number R* m — 15.63), x y = (Petschek-like solution). As R* m increases, x y also 
increases. In the limit of large R* em , x y approaches to unity. 
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